Theory of a one-dimensional double-X atom interferometer 
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The dynamics of an atom waveguide X-junction beam splitter becomes truly ID in a regime of 
low temperatures and densities and large positive scattering lengths where the transverse mode 
becomes frozen and the many-body Schrodinger dynamics becomes exactly soluble via a generalized 
Fermi-Bose mapping theorem. We analyze the interferometric response of a double-X interferometer 
of this type due to potential differences between the interferometer arms. 
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I. INTRODUCTION 



Recent advances in atom de Broglie waveguide tech- 
nology jl], |2[ ||, @j Sand its potential applicability to 
atom interferometry |(| and integrated atom optics [[| 0| 
create a need for accurate theoretical modelling of such 
systems in the low temperature, tight waveguide regime 
where transverse excitations are frozen out and the quan- 
tum dynamics becomes essentially one-dimensional (ID) 
(Tonks-gas limit). It has been shown by Olshanii ||, 
and also by Petrov et al. ||, that at sufficiently low 
temperatures and densities and high transverse frequen- 
cies loq, where thermal and longitudinal zero-point en- 
ergies are small compared with Huq, the transverse de- 
grees of freedom of a Bose-Einstein condensate (BEC) in 
an atom waveguide are frozen in the ground transverse 
mode and the dynamics becomes one-dimensional. If, in 
addition, the scattering length is large and positive, the 
dynamics reduce to those of a ID gas of hard core, or 
impenetrable, point bosons ||, This is a model for 
which the exact many-body energy eigensolutions were 
found in 1960 using an exact mapping from the Hilbert 
space of energy eigenstates of an ideal gas of spinless 
fermions to that of many-body eigenstates of hard core, 
and therefore strongly interacting, bosons Jl(], Re- 
cently two of us have extended this mapping technique in 
order to treat the exact many-body dynamics of Tonks 
gases p2| , |l3| and the exact many-body ground-state of 
a trapped Tonks gas Q . 

"Waveguide on a chip" technology (l^, [l6| has already 
advanced to achievement of BEC on a magnetic surface 
microtrap and construction of a Y-beam splitter on a 
chip Jl8| , and the theory of a multimode double-Y inter- 
ferometer has been developed Jl9[ ] in an approximation 
where atomic interactions are neglected. Here we are in- 
terested in the the opposite limit of the Tonks-gas regime 
[^|, ^| where interatomic interactions and two-body cor- 
relations play a major role. Such a study is motivated 
by our previous demonstration jil| that fermionization 



strongly inhibits interference in an adiabatically split and 
recombined Tonks gas. The Fermi-Bose mapping and dy- 
namics of a pulsed interferometer arc more complicated 
and it is not clear a priori whether significant interference 
effects will occur in this case. This motivates the present 
investigation. 

We set up our model for the double-X interf erome ter 
in Section II as an effective ID problem. 
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In Sec. 

develop a generalized mapping theorem and thereby map 
the T onks state to a ID Fermi gas of free fermions, in Sec. 
IIIB we define the boundary conditions at the junctions 
in terms of single particle wave functions, and in the re- 
mainder of Sec. |H we examine the limits of validity of 
the model and we develop a general expression for the 
bosonic spatial density profiles in terms of simpler Fermi 
orbitals which are not subject to the X-splitter boundary 
conditions. Analytical expressions for the interferometric 
response are developed in Sec. fv|and numerical results 
are presented in Sec. [v|. 



II. DOUBLE-X INTERFEROMETER 

The interferometer model of Andersson et al. |l9| ] 
consists of an atom waveguide which is first split into 
two identical waveguides by a symmetrical Y-junction, 
followed by a reverse Y-junction which recombines the 
output of the two arms into a single waveguide, the 
"exit port" , where interference fringes appear as a re- 
sult of potential differences between the two interferome- 
ter arms. The emphasis there was on effects of multi- 
transverse-mode propagation, assuming the transverse 
trap frequency low enough and/or longitudinal density 
high enough and/or temperature high enough that many 
transverse modes of the guides are energetically accessi- 
ble and many-body correlations are negligible. In the op- 
posite Tonks limit where the transverse wave function is 
frozen, the dynamics fermionizes and multiple longitudi- 
nal modes play a crucial role, due to the Fermi-Bose map- 
ping theorem |u| [H], [l^, [l3| |l4j which introduces a Fermi 
sea of many longitudinal orbitals into the many-boson 
dynamics. Rather than the double-Y geometry of An- 
dersson et al, we assume a double-X geometry since its 
greater symmetry and reversibility makes implementa- 
tion of unitary straightforward; It is essentially the same 
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FIG. 1: Schematic diagram of the interferometer. The inset 
shows that the junctions are actually avoided crossings with 
tunnelling between the waveguides. 



model used by Scully and Dowling [^0[ in their treat- 
ment of an atom waveguide analog of a Mach-Zehnder op- 
tical interferometer. We assume that the entrance arms, 
upper and lower arms after the X-splitter, and exit arms 
after the X-recombiner all satisfy the necessary and suf- 
ficient conditions for the Tonks regime ^ [)) . We can re- 
tain an effective ID description despite the splitting and 
recombination (a) by assuming that corresponding up- 
per and lower interferometer arms have the same length, 
so that the same longitudinal coordinate x can be used 
for both arms and (b) by introducing a pseudospin rep- 
resentation with s = 1 (— 1) labelling atoms in the upper 
(lower) interferometer arm. Our basic model of a double- 
X interferometer is shown in Fig. ^, and we remark that 
the same model applies to an atomic loop mirror em- 
ploying an X-input-coupler. The inset indicates that the 
two apparent X-crossings are really "avoided crossings" 
of the waveguides, with tunneling between the two in 
their regions of near tangency. We idealize these tun- 
neling regions to points, see discussion following Eq. (||), 
taking the X-splitter to be at x — L\ and X-recombiner at 
x = L 2 and treating the effects of tunneling by boundary 
conditions specifying the wave function discontinuities at 
these two points. Then x specifies the position of an atom 
in one of the entrance arms if x < L\, one of the inter- 
ferometer arms if L\ < x < L 2 , and one of the exit arms 
if x > L 2 . Details of these boundary conditions will be 
given in the following section. 



III. THEORETICAL MODEL 

A. Generalized Mapping Theorem 

Our objective is to find dynamical solutions for an 
TV-boson Tonks gas in the X-splitter geometry which 
are exact solutions of the time-dependent many-body 
Schrodinger equation (TDMBSE) except when any of the 
N particle positions Xj is equal to L\ or L 2 (X-junctions), 



satisfy a point hard core impenetrability constraint, have 
Bose symmetry (totally symmetric), and include the 
effect of tunneling at the X-junctions via appropriate 
boundary conditions. We do this by generalizing our pre- 
vious Fermi-Bose mapping theorem 

El His to 

the present more complicated geometry. 

A system of N atoms is described by solutions 
y(x 1 ,s 1 ;---;x N ,s N ;t) of the TDMBSE = ihd^/dt. 
In the Tonks limit (large scattering length, low linear den- 
sity, tight transverse confinement) ||, ||] the two-particle 
interaction behaves as a hard core of vanishing diameter 
a — ► 0. This is conveniently treated as a constraint on 
allowed wave functions: 

* = if Sj = s k and x 3 ■ = x k , 1 < j < k < N, (1) 

implying that atoms in the same arm cannot inter- 
penetrate whereas those in different arms have no 
such constraint. We start from fermionic solutions 
^f(xi, si; ■ ■ ■ ; xn, sn; t) of the TDMBSE which are an- 
tisymmetric under all space-pseudospin pair exchanges 
(xj,Sj) <-> (x k ,s k ), hence all permutations. Generalizing 
the previous definition (ilj, |l| |l|, || Q, define a "unit 
antisymmetric function" A by 

A(x 1 ,s 1 ; ■ ■ ■ ;x N ,s N ) = a(xj, Sj\ x k , s k ) 

l<j<k<N 

a(Xj,Sj;x k ,Sk) = 5 Sk:Sj sgn(x k - x } ) 

+ SsjiSsk-i - $ Sj -i$ Sk i, (2) 

where sgn(x) = +1(— 1) if x > Q(x < 0). For a given 
antisymmetric ^f, define a bosonic wave function 'Fa 
by 

^b(xi,si; ■■■ ;x N ,s N ;t) = A(xi,si; ■ ■ ■ ;x N ,s N ) 

x ty F (xi,si;---;xff,a N ;t). (3) 

It satisfies the impenetrability constraint (|l|) because any 
spinlcss fcrmion wave function is space-antisymmetric 
under exchange of like-pseudospin pairs. It is totally 
symmetric (bosonic) under space-pseudospin permuta- 
tions, obeys the same boundary conditions as and 
the Bose TDMBSE H^ B = ihdH B /dt follows from the 
Fermi one H^p = ihd^F/dt by staightforward gener- 
alization of the previous argument O, [l^, |l^, fl4| , 
under the assumption that the Hamiltonian is the sum 
of kinetic energy operators and ordinary (non-operator) 
potentials. These potentials may depend on pseudospin, 
i.e., they may differ between the left and right interfer- 
ometer arms, as is the case in modelling interferometric 
detectors. Since the only interatomic interaction in the 
Tonks limit is the point hard core treated by the con- 
straint (0), and this constraint is satisfied automatically 
by the iV-fermion wave functions ^f in @ as a con- 
sequence of their antisymmetry, it follows that $f is a 
dynamical ideal Fermi gas state which can be written as 
a Slater determinant 



^f(xi, si; • • • ;x N , s N ;t) 



1 N 

det^ cj) n (xj,Sj;t), 



/JVT n,j=l 



(4) 
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where <p n are N orthonormal solutions of the single- 
particle TDSE. 



B. X-splitter boundary conditions 



The TDMBSE is 



h 2 d 2 

2m dx 2 



+ V s (xj,t) -m 



* = 0, (5) 



where V s (x, t) are longitudinal potentials assumed to be 
nonzero only within the interferometer arms L\ < x < 
L2. The mapping theorem guarantees that ^b will sat- 
isfy this TDMBSE at all points (x\, s±; ■ ■ ■ ; xn, sn) where 
^>f does. However, and hence ^b fail to satisfy the 
TDMBSE at the points Xj = L\ and Xj — L2 where tun- 
neling introduces discontinuities in the wave function; at 
these points the TDMBSE is replaced by boundary con- 
ditions specifying the discontinuities. Each of the orbitals 
in the fermionic solution (|) of the TDMBSE (|) satisfies 
the single-particle TDSE 



h 2 d 2 „, . B 



<p„{x, s;t) = 0, (6) 



everywhere except at the X-junctions x — L\ and x = L2, 
where this TDSE is replaced by a boundary condition. 
Since evolution according to the TDSE is unitary, the 
transformation from the left to the right side of each 
junction can be implemented by a unitary 2x2 matrix. 
Assume that at the junction half of the probability ini- 
tially in the upper arm remains in the upper (s = 1) 
arm and half tunnels to the lower arm, and similarly for 
probability entering from the lower arm, and introduce a 
Pauli pseudospinor notation 



<f) n (x,+l;t) 
(f> n (x,-l;t) 



(7) 



Then a simple boundary condition which implements the 
splitting at the first X-junction is pO, p2[ 



* n (Li+,t) =T*„(L 1 -,t) ! 



(8) 



where L\— (£1+) denotes the left (right) side of the junc- 
tion and T is the unitary matrix 



V2 



1 -1 
1 1 



(9) 



Thus the amplitudes which remain in the upper or lower 
arms as they pass through the avoided crossing junction 
suffer no phase shifts, that which tunnels from the upper 
to the lower arm suffers a phase shift of 7T, and that 
which tunnels from the lower to the upper arm suffers 
no phase change. It can be shown that this matrix can 
be written in the form e lV — (1/V2)(1— o++cr_) in terms 
of an interaction V = iA(<7+ — <r_) where &± are Pauli 



pseudospin raising and lowering operators implementing 
the tunneling and A is taken equal to 7r/4 in order that 
probability splits equally between the upper and lower 
arms. More general phase shifts could be used subject 
to the proviso that T must be unitary, but the simple 
choice here is sufficient. At the recombiner (x = L2) the 
same matrix boundary condition is applied once more to 
generate the amplitudes in the upper and lower output 



C. Limits of validity 

In our simple model the tunneling regions giving rise 
to the X-splitters are replaced by point boundary con- 
ditions, and we now examine the conditions for this to 
be valid. In particular, in the boundary condition (||) we 
treat each orbital separately, meaning that the tunneling 
is assumed to proceed as for single particles and to not 
be modified by many-body interactions. For a 50:50 X- 
splitter, and assuming a tunneling region length L coup , 
the energy splitting AE << Hloq between the symmet- 
ric and anti-symmetric transverse modes of the tunnel- 
coupled waveguides must obey (AE/h) ■ (L coup /v x ) = 
7r/4, where v x is the longitudinal velocity of the atoms 
entering the X-splitter. Here L coup / v x is the time of flight 
of the atoms through the tunneling regions, and the 7r/4 
phase-diffcrcncc between the modes is that required for a 
50:50 splitter. For single-particle tunneling to be valid we 
require that the energy splitting is large compared to the 
many-body energy per particle of the incoming Tonks gas 
of mean density p : AE » h 2 n 2 p 2 /6m Q |f). When 
this condition is satisfied the atoms tunnel individually 
and the transfer characteristics of the X-splitter should 
therefore be the same for bosons and fcrmions. For an ini- 
tial gas in its ground-state which is harmonically trapped 
with longitudinal frequency to we have p ~ N /xf, where 
xp = \/2Nx osc , x osc = \J h/muj being the single particle 
ground-state width ^3|. In order that the tunneling re- 
gion not create large density gradients in the atomic gas 
we also require L coup > xf , that is the coupling region is 
longer than the width of the initial trapped gas. Further- 
more, under the assumption L coup << {L2 — L\) we may 
treat the tunnel regions as points. Combining these re- 
sults together, we obtain the following conditions for the 
action of the X-splitters to be treated by the boundary 
condition (||) 

v x » N ■ 

-^coup 

■ijj (10) 
(L2-L1) > L coup > V2Nx osc . (11) 

We shall hereafter assume these conditions arc satisfied. 



D. Construction of interferometer solutions 

The boundary conditions (||) at the X-splitters can be 
implemented as follows: Let u n (x,s;t) be solutions on 
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x e [—00,00] of the same sing le-particle TDSE © as 
the tf> n (x,s;t) satisfy, but now no boundary condition is 
imposed on these orbitals, so they are continuous with 
continuous gradient at all x including x — L\ and x = L 2 . 
Instead, an initial condition 



u n (x, s;t = 0) = u n {x), 



(12) 



is imposed. Note that this initial condition is indepen- 
dent of s (the same in both entrance arms). Here t = 
is a time shortly after loading of the N atoms into the 
entrance arm and the u n (x), defined for 1 < n < N, 
are a set of orthonormal orbitals which are negligible at 
t = for x > L\. Then for t > and x < L%, define new 
orbitals </>„ by 



cj) n (x,s;t) = u n (x,s;t)5. 



si 



x < L 



1 • 



(13) 



These </>„, which in our model represent the actual or- 
bitals occupied by the atoms, satisfy the same TDSE (|^) 
as the u n for x < L\, but they are nonvanishing only in 
the upper arm [s = 1), corresponding to our assumption 
that the atoms are to be loaded only into the upper arm. 
On the other hand, the u n are nonzero in both the upper 
and lower arms as a consequence of their initial condi- 
tions (||) and TDSE (§). This allows us to first solve 
the TDSE without imposing boundary conditions at the 
junctions x — L\ and x = L%, then construct solutions 
(f) n of the TDSE satisfying the proper boundary condi- 
tions by defining them in terms of the u n by definitions 
which are properly discontinuous at the junctions. The 
effect of the boundary condition at x — L\ is such that 
the upper-arm amplitude u n (L\, l;t) splits equally, with 
prefactors l/v2, between the upper and lower arms at 
the junction. Note that u n {Li, 1; t) =u n {L\,— l;t) since 
V s vanishes for x < L\ and the initial condition (|l2|) is 
independent of s. It then follows that the solutions of 
the TDSE in the interferometer arms are simply 



*n(ar,*) = -7= 



u n (x, 1; t) 
u n (x,-l;t) 



Li < x < L 2 (14) 



In general these u n depend on s = ±1 since the potentials 
V s in the arms L\ < x < L 2 do. 

At the recombiner x — L2 the splitting matrix T is 
applied once more, resulting in the solutions 



u n (x,l;t) - u n (x,-l;t) 
u n (x,l;t) +u n (x,-l;t) 



x > L 2 (15) 



in the output arms, where the potentials V s are again 
zero. 

The above equations uniquely and straightforwardly 
determine the many-fermion wave function typ in terms 
of N orbitals <j> n . Then within the limits of validity of 
our theory given in Sec. ( Ill C] ), under which the tunnel- 
ing occurs as for single-particles and is hence the same 
for bosons and fermions, the many-boson wave function 
is given by \fs = Atyp. Since A 2 — 1, this has the 
important consequence that the single-particle densities 



p(x, s;t) are the same for the Bose Tonks gas as for the 
ideal "spinless" Fermi gas of the mapping theorem, al- 
lowing us to compute Bosc dynamics and interference 
fringes for this model of a Tonks gas interferometer in 
terms of the corresponding solutions for the ideal Fermi 
gas. This will be done in the following sections. In addi- 
tion to providing the basis for an exact calculation of the 
behavior of a Tonks gas of bosonic atoms in this interfer- 
ometer model via the Fermi-Bose mapping theorem, this 
fcrmionic calculation is of interest in its own right, since 
a spin-polarized Tonks gas of fcrmionic atoms is dynami- 
cally equivalent to the same ideal Fermi gas inasmuch as 
it satisfies the impenetrability constraint (Q) automati- 
cally as a consequence of spatial antisymmetry. 



E. Construction of spatial density profiles 

It follows from orthonormality of the u n [x) 
and the unitarity of Schrodinger evolution that 
the u n (x, s;t) satisfy the orthonormality relation 
w* (x, s; t)u m (x, s; i)dx = 8 nm for each fixed s = ±1. 
One can then show with Eqs. ( |ll|]l5| ) that the spinor 
orbitals 3> ra are orthonormal in the usual sense 

& n (x,t)& m (x,t)dx 

-00 

/OC 
4>* n (x,s\t)<j) m (x,s;t)dx 

1 f 00 

= - ^ / u* n {x, s;t)u m (x,s;t)dx = 8 nm , (16) 

in spite of their discontinuities at the X-junctions x = L\ 
and x = L 2 . Note that the orthonormality of the u n 
holds for each fixed s (separately in each arm), whereas 
the orthonormality relation for the 3>„ involves both sum- 
mation over s and integration over x. This is to be ex- 
pected physically, since the X-junctions mix probability 
density between both arms. Using the Slater determinan- 
tal expression for ^p in terms of the 4> n (x,s;t) and the 
fact that l^sl 2 = I^fI 2 , it then follows that the single- 
particle density for both the Bose and Fermi systems is 



JV-l 



p(x,s;t) = 22 \<f>n(x,s;t)f 



(17) 



To recapitulate, to generate the density profiles we first 
solve the TDSE for each u n (x, s,t),n — 1, . . . , N in both 
arms, which does not involve the X-splitter boundary 
conditions. Next the orbitals (f> n (x,s;t) at the detector 
are constructed from the boundary condition JIe|) , and fi- 
nally the density profile is obtained from Eq. (|17|) above. 



IV. ANALYTICAL SOLUTION 

Hansel et al. Jl6|] have demonstrated a magnetic con- 
veyor belt for atoms that could be used to launch an 
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atomic gas into interferometer, and recently transport 
of BECs by optical tweezers has been achieved by Gus- 
tavson et al. |24|], who suggested that this technique 
could be used to load a BEC into a "waveguide on a 
chip" . We simulate the tweezer loading method by choos- 
ing the iV-fermion wave function in the entrance or upper 
arm arm s = +1 in Fig. [l] as the ground-state of a har- 
monic trap of frequency u> with an overall translational 
velocity v x = hk/m, k being the longitudinal wavevector. 
The time when the Fermi wavepacket is released from the 
optical trap into the entrance arm is taken as the time 
origin t = and the position of the packet-center at that 
instant is taken as the space origin x — 0. Then, the or- 
thonormal orbitals defining the initial condition in fll2] ) 
are conveniently chosen to coincide with the lowest ./V 
Hermite-Gaussians 



t (x) = C n e~ x ^ 2x ^H n 



ikx 



C n — 



1 



(18) 



where the N-fermion ground-state has one particle in 
each of the first N modes, and x osc = yjK/moj as be- 
fore. The Fermi- wavev ector associated with this N-atom 
ground-state is k F = y/2N /x osc j25|, [2(|, and this serves 
as a measure of the wavevector on the highest occupied 
plane- wave orbital. We assume throughout that k » kp 
so that the center-of-mass velocity v x of the atomic 
wavepacket exceeds the expansion of the wavepacket af- 
ter it is released. 

Our goal in this paper is to examine the response of 
the double-X interferometer to a variety of perturbations 
V s ■ An interesting choice of perturbing potentials in the 
interferometer arms is one which is a delta-function pulse 
in time, V s (x,t) = HA s (x)S(t—to)- Then, the orbitals just 
after such a perturbation are related to those just before 
by a spatially varying phase-shift A s (x): 



(19) 



l (x,s;t ) , L l < x < L 2 

Such an approach was used by Rojo, Cohen, and Berman 
to study the response of a Tonks gas to a delta-pulsed 
optical lattice potential p{j| , and we have used it more 
recently to investigate gray solitons in a Tonks gas 

El- it 

is an idealization to zero pulse length of phase imprinting 
techniqes used in a number of recent experiments [^7| . |2q . 

As a simple but concrete example, suppose that an 
external force acting on the interferometer exerts a dif- 
ferent impulsive force F s = sh[k p + 2Qx]S(t — to) on each 
arm at a time to at which the wavepackets are located 
within the interferometer. (We split the impulsive force 
between the two arms for ease in the analytic solution 
but the same physics arises if the force is applied to one 
arm). The force term proportional to k p represents a spa- 
tially uniform force, whereas the term proportional to Q 
represents a spatially inhomogeneous force. The corre- 
sponding potential is V s (x,t) = HA s (x)5(t — to) with 



where ip is a phase difference between the two interfer- 
ometer arms due, for example, to the Sagnac effect if 
the interferometer is rotating. The factor s = ±1 means 
that the impulse is oppositely directed in the left and 
right arms. In the absence of a quadratic phase varia- 
tion, Q — 0, via the equivalence principle this can be 
viewed as a simple model of a detector of rotational ac- 
celeration, an "atom gyroscope" . 

The details of the time evolution of each orbital 
u n (x,s;t) in either arm from the time of loading un- 
til the time of detection are given in the Appendix: 
The expressions for these orbitals are somewhat cum- 
bersome, but in general they can be written in the form 
u n (x,s;t) — U n (x, t) exp(i^„(x, s, t)), where U n are real 
Hermite- Gaussian functions and £ n are phase factors, so 
\u n \ 2 = \U n \ 2 - The time-evolution is the same in each 
arm until the time t = to when the pulses are applied. 
That time is chosen so that the pulses are well con- 
tained in the interferometer arms and the amplitude of 
the wavepacket at the junctions is neglig ible (see Fig. |). 
Thus for subsequent propagation the extra phase ( |l9| ) ac- 
quired by each orbital due to the pulsed potentials can 
be taken to extend over all values of x although the po- 
tentials are limited to the lengths of the interferometer 
arms only. 

If the potentials have no quadratic spatial dependence 
Q = 0, the phase acquired is the same for each orbital, 
as seen from Eq. (A5) in the Appendix. In that case 
the densities in the upper arm (s = +1) and the lower 
arm (s = —1) can be expressed in terms of real- valued 
orbitals U n centered at x± = xq + h(k ± k p )(t — to)/m, 
see Eq. 

1 N-l 

p(x,±l;t) = ^J2 [Un(x-x+,t) + U 2 (x-x_,t)] 



1 



JV-1 
\ ^ 



T-cos[0(z,t)] ^ U n {x - x+,t)U n {x - x.,t){2\) 

n=0 

where xq = v x to is the position of the center of the pack- 
ets at the scaled time of application of the pulse tq = uto, 
and x± the centers of the wavepackets in each arm at the 
time of detection r = Lot. The spatial phase-modulation 
9 due to interference and the period A of that modulation 
are given by 



(x, t) = ifi + 2k p 



l+r 2 

1 + TQT ' 



1+TQT 
1+T 2 



x 



(22) 



A s (x, t) = -s[<p/2 + k p x + Qx 2 ]5(t - t ), 



(20) 



where x m = {x + + xJ)j2 is the mean position of the 
recombined wavepackets at the time of detection. The 
exact analytic solution ( pT| ) will be used in the following 
sections to evaluate the response of the double-X inter- 
ferometer to the input atom wavepacket under varying 
applied perturbations. 

Equation ( pl| ) for the density profile in each arm has an 
illustrative limiting case. For small applied wavevectors 
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k p <C k, the centers of the wavepackets do not separate 
much x± w v x t, and we have U n (x — x+,t) ~ f/„(a; — 
:r_ , i) . In this limit the density in each arm becomes 



p(x,±l;t) ~ - [1 Tcos[6»(x,] 



AT-1 

E 

n=0 



0*(a!-a; + ,i), (23) 



this expression being exact for fc p = 0. The factor 
^2 Un(x — x+, t) is the density profile of the expanding 
Tonks gas, and is identical in the present approxima- 
tion to that in the absence of the applied perturbation. 
Thus, for k p <gi k the response of the interferometer to 
the perturbation is in the spatially dependent phase fac- 
tor 8(x,t). For the case k p = we have 0(x,i) = tp, 
so the densities in each arm do not show spatial fringes, 
and we may integrate the density over x to find exact 
expressions for the number of atoms N± in each output 
arm 

N 

N± = - [1 T cob(vj)] . (24) 

Then if we set tp = tpb+Atp, where ipb is a bias phase-shift, 
created, for example, using an optical dipole potential 
applied to one arm, and Atp is a phase difference between 
the two arms due to the Sagnac effect, we find (N + — 
N-) = — N cos(ipb + Atp). Setting tpb = it/2 we obtain 



(JV+ -N-) = Nsm(Aip). 



(25) 



Thus, measuring the difference in numbers of atoms ex- 
iting the two interferometer arms provides a measure of 
the excess phase-shift Aip = 2mujjiA/h due to rotation, 
for example, where lor is the rotation frequency and A 
the enclosed area of the interferometer. 

When fcp^O Eq. ( |2l| ) shows that the output density in 
each arm can exhibit spatial fringes with period A as the 
phase 6{x, t) becomes x-dependent. The fringes become 
apparent when 2xf/\ = 2^/2Nx scl 'A > 1, so that there 
is more than one fringe under the atomic wavepacket 
width. The condition on the perturbing wavevector to 
observe fringes is therefore 



> 



1 



1 



) \l+T T 



(26) 



However, for momentum kicks k p such that the distance 
(x + — X-) = 2hk p (t — to)/m between the wavepacket 
centers in the two arms becomes greater than the width 
Xpyfl + t 2 of the expanding Tonks gas at scaled detec- 
tion time t = ujt, then the fringes will tend to vanish 
as the exiting wavepackets from the two arms no longer 
overlap. This yields the following condition for the exit- 
ing wavepackets to not overlap 



k 



> 



N X osc yJl + ■ 



(27) 



2 (x m - x ) 

where x m is the mean position of the recombincd 
wavepackets at the scaled time r = cut of detection. 
Furthermore, if r ?a Lu(x m /v x ) <C 1, then the atomic 
wavepacket will not spread much during its passage 
through the interferometer. 
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FIG. 2: The shape of the wavepacket in either arm before the 
pulses are applied, shown when the packet center is at a;o = 2 
mm . 



V. NUMERICAL RESULTS 

In this section we present numerical solutions from our 
model to illustrate the interferometric response of the 
double-X atom interferometer to a variety of imposed 
perturbations. These illustrations are intended to display 
some interesting features of the theory and not to serve 
as detailed study of a particular mode of inteferometer 
operation, eg. a rotation sensor. 

For our numerical study of the interferometric re- 
sponse of the double-X interferometer we choose param- 
eters corresponding to the optical-tweezer loading ex- 
periment of Ref. p4j , with longitudinal trap frequency 
v = 4Hz => u) = 25.1 Hz, which gives an oscillator length 
of x osc = 1.04 x 10 -2 mm for sodium atoms. The optimal 
transfer velocity of 70 mm/s is chosen to be the initial 
velocity v x of the wavepacket launched into the interfer- 
ometer, which corresponds to a wavevector k = 267x~} c . 
The arm-lengths in the interferometer are fixed by tak- 
ing L\ = 0.5 mm, L2 — 2.5 mm and the region of de- 
tection to be centered at Xd — 3 mm, 0.5 mm beyond 
the X-recombiner, see Fig. |[ Figure || shows the den- 
sity in each of the two interferometer arms for N = 51 
atoms just before the perturbing pulse is applied, and the 
wavepacket is seen to be well contained between the two 
X-splitters at L\ and L2. In the following discussion the 
output atom densities are calculated around the point of 
detection x m = 3 mm. 

The response of the atom interferometer to a net im- 
posed phase-shift <p between the two arms has already 
been addressed in the discussion leading to Eq. (p5|) , so 
here we set tp = and look at the effects of an applied 
force. Figure § shows the output density profiles in the 
upper (s = +1) and the lower (s = —1) arms of the in- 
terferometer for an initial wavepacket of N = 51 atoms 
moving with velocity v x — 70 mm/s, pulsed force ap- 
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2.5 3 3.5 

x (mm) 



FIG. 3: Density profile of the recombined wavepackets in the 
upper (s = +1) and the lower (s = —1) arms for an initial 
wavepacket of 51 atoms moving with velocity v x = 70 mm/s, 
and the pulsed potentials V s (x,t) applied when xq — 2 mm: 
(a) kp/k = 10" 4 and Q = 0, (b)k p /k = 10~ 3 and Q = and 
(c)k p /k = 10" 3 and Q/k = 5 x 10" 4 mm" 1 . 



plied at xo = 2 mm with (a) k p /k = 10 4 and Q = 0, 
(b)h p /k = 10~ 3 and Q = and (c)fcp/fc = 10" 3 and 
Q/k = 5 x 10~ 4 mm" 1 . For a spatially uniform force 
(Q = 0) the condition (|2q ) to observe fringes yields 
fcp/fc > 0.6 x 10 -3 , and this is consistent with case (a) 
which shows no fringes k p = 10~ 4 , and case (b) for 
k p = 10~ 3 which shows a couple of fringes. Thus, the 
period of the density fringes in each arm can act as a 
measure of the wavevector of the applied force via Eq. 
(p2|). However, as shown in Fig. |3|(c), which is the same 
as (b) except with Q ^ 0, if the force has a spatially 
non-uniform component the visibility fringe can be re- 
duced. We have found that the density fringes are sensi- 
tive to small spatially inhomogeneous compenents of the 
applied force. As the wavevector k p of the applied per- 
turbation is further increased beyond that used in Fig. 
^(b), again with Q = 0, the number of density fringes 
increases but eventually their contrast decreases. In par- 
ticular, according to the condition (p7j) if k p /k > 0.05 
then the exiting atom wavepackets no longer overlap and 
interference fringes are not possible. 

As a second example, Fig. H shows the response of 




FIG. 4: Number of atoms in the upper (s = +1) and the 
lower arms as a function of linear pulse strength k p , with 
Q — and with the time of application of the pulses V s (x,t) 
fixed at when the packet centers are at xo = 2 mm. The 
initial wavepacket has 51 atoms moving with velocity v x — 70 
mm/s. 



the interferometer to a force of varying wavevector k p for 
Xq = 2 mm and N = 51 as before. (Keeping the position 
or time of application of the applied force fixed is some- 
what artificial but serves to illustrate our theory). In par- 
ticular, the figure shows the number of atoms exiting each 
arm as a function of k p . For k p /k -C 0.6 x 10 -3 , so that 
there are no density fringes under the atom wavepacket, 
the atom numbers exhibit almost complete periodic oscil- 
lations. We can understand this behavior from Eq. ( p3) 
for the densities in each arm: When k p /k <C 0.6 x 1CP^ 
the phase 0(x,t) varies little over the spatial extent of 
the atom wavepacket, in which case the integrated atom 
numbers in each arm are given to a good approxima- 
tion by Eq. (|24|) with ip replaced by 2k p xo, which shows 
periodic exchange. In contrast, for k p /k > 0.6 x 10~ 3 
the atom number oscillations are much smaller in ampli- 
tude and centered around 50:50 splitting. In the regime 
k p /k ^S> 0.6 x 10 -3 the phase 9(x,t) varies significantly 
over the spatial extent of the atom wavepacket, and upon 
integrating Eq. ( |2^ ) over x to obtain the atom numbers 
in each arm, the cosine terms will tend to integrate to 
zero due to their oscillatory nature, giving 50:50 split- 
ting. 



VI. SUMMARY AND CONCLUSIONS 

In summary, by introducing a generalization of the 
Fermi-Bose mapping for coupled waveguides we have de- 
veloped the theory of a double-X atom interferometer in 
the Tonks-gas limit ^] of tight transverse confinement 
and impenetrable interactions between the atoms. We 
believe this is a significant development in view of cur- 
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rent efforts to create ID integrated atom interferometers, 
and in waveguide on a chip technology in general, for iner- 
tial and other sensors. In particular, we have calculated 
the interferometric response of the double-X device for 
a variety of imposed external perturbations to illustrate 
the utility of our theory. A key result of our work is to 
show that despite the fact that the we are in the Tonks- 
gas limit, where the condensed fraction of the initial gas 
tends to zero for large number of atoms, the interferome- 
ter can still show high visibility density fringes at the out- 
put. That is, our results show the ID waveguide devices 
with ultracold but non-condensate atom sources can still 
exhibit a high degree of first-order coherence. Further- 
more, due to the Fermi-Bose mapping our results apply 
equally well to a spin-polarized fermionic atom source, 
for which s-wave scattering is absent. 
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APPENDIX A: TIME EVOLUTION OF THE 
ORBITALS 



where the new center of the wavepacket is at xo = kto- 
After the pulsed potentials V s (x,t) are applied the or- 
bitals acquire an additional space-dependent phase dif- 
ferent in each arm 

u n (x, s; t+) = e is ^+^^u n (x, s; t„ ). (A3) 
Evaluating the subsequent free propagation till a time 
t > to after the wavepacket from each arm has completely 
passed through the recombiner, one finds the functional 
form of the modes coming from the two arms to be 



u n (x,±l;t) = 



C n e—^ 



^(i±2nT d ) + i(T±2 To nT d ) 
(i±2nT d ) + i{T±2T nT d )\ 

- [Xo + («± ± 2 Xo ^)T d ] \ 



V[l±2fir ti ] 2 + [r±2Tor!r d ] 2 



(A4) 



with k± = k± n p , and the argument of the exponential 
is given by 

F±(x,t,t ) = -(l±2fh- d )xo - n ± T d [2 Xo + in±(l + ir a )} 
+2x[xo + i«±(l + no)] - [(1 ± 2f>T ) - i(±2n)} X 2 
G±(X, r, r ) - 2[(1 + it) + (±2fi)r d (l + ir )]. 



In the absence of external potentials the time evolved 
orbitals u n (x, s; t') at time t — t' + At are obtained from 
those at time t by applying the retarded free particle 
Green's function 



u n (x,s;t) = 



V2TTiujAt 



f 

J — ( 



dx' u n (x',s;t')e '^Ja* (Al) 



For brevity in notation notation we introduce the scaled 
variables To — toto for the time of application of the pulse, 
t = ujt for the time of detection, and T d = uj(t — to)- 
We also introduce Greek letters to designate the scaled 
variables x = x/x osc and k = kx osc and £1 = Qx osc . The 
free propagation of each orbital u n (x, s; t) via either arm 
is identical up until the instant t = t r before the pulses 
are applied 



u n (x, s; to ) 




-(x-xo) 

xe 1KX e 2 d+-o) H, 



(A2) 



When the potential has only a linear dependence on x 
i.e. — this expression acquires a more transparent 
form 



u n (x,±l;t) 



(K ± T d + K To) 



1 — IT 



l + ir 



n/2 



" °" - (x-(XQ + ^r d y (A5) 

VI +T 2 



2(l+iT) 



This expression contains a phase-independent factor 
which is a Hc rmitc-G aussian centered at x± = x a + K±T d 
with width Vl + t 2 , that we represent in terms of un- 
sealed variables 



-i*-*±i 2 

U n (x-x±,t)=— — , .,„,, - H, 



x — x± 



[1+M) 2 ]V4 »L 



. (A6) 
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